home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slantp.z / slantp
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAANNNNTTTTPPPP((((3333FFFF))))                                                          SSSSLLLLAAAANNNNTTTTPPPP((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLANTP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a triangular
  11.      matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      REAL FUNCTION SLANTP( NORM, UPLO, DIAG, N, AP, WORK )
  15.  
  16.          CHARACTER DIAG, NORM, UPLO
  17.  
  18.          INTEGER   N
  19.  
  20.          REAL      AP( * ), WORK( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      SLANTP  returns the value of the one norm,  or the Frobenius norm, or the
  24.      infinity norm,  or the  element of  largest absolute value  of a
  25.      triangular matrix A, supplied in packed form.
  26.  
  27.  
  28. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  29.      SLANTP returns the value
  30.  
  31.         SLANTP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  32.                  (
  33.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  34.                  (
  35.                  ( normI(A),         NORM = 'I' or 'i'
  36.                  (
  37.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  38.  
  39.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  40.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  41.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  42.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  43.  
  44.  
  45. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  46.      NORM    (input) CHARACTER*1
  47.              Specifies the value to be returned in SLANTP as described above.
  48.  
  49.      UPLO    (input) CHARACTER*1
  50.              Specifies whether the matrix A is upper or lower triangular.  =
  51.              'U':  Upper triangular
  52.              = 'L':  Lower triangular
  53.  
  54.      DIAG    (input) CHARACTER*1
  55.              Specifies whether or not the matrix A is unit triangular.  = 'N':
  56.              Non-unit triangular
  57.              = 'U':  Unit triangular
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAANNNNTTTTPPPP((((3333FFFF))))                                                          SSSSLLLLAAAANNNNTTTTPPPP((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      N       (input) INTEGER
  75.              The order of the matrix A.  N >= 0.  When N = 0, SLANTP is set to
  76.              zero.
  77.  
  78.      AP      (input) REAL array, dimension (N*(N+1)/2)
  79.              The upper or lower triangular matrix A, packed columnwise in a
  80.              linear array.  The j-th column of A is stored in the array AP as
  81.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  82.              if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.  Note
  83.              that when DIAG = 'U', the elements of the array AP corresponding
  84.              to the diagonal elements of the matrix A are not referenced, but
  85.              are assumed to be one.
  86.  
  87.      WORK    (workspace) REAL array, dimension (LWORK),
  88.              where LWORK >= N when NORM = 'I'; otherwise, WORK is not
  89.              referenced.
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.